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Abstract 

In this work, a model to study the coupUng between a semiconductor qubit and two time- 
dependent electric fields is developed. By using it in the resonantly monochromatic double dressing 
regime, control of the local density of optical states is theoretically and numerically demonstrated 
for a strongly confined exciton. 

Drastic changes in the allowed energy transitions yielding tunable broadening of the optically 
active frequency ranges, are observed in the simulated emission spectra. The presented results are 
in excellent qualitative and quantitative agreement with recent experimental observations. 
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1. INTRODUCTION 



Few decades ago, artificial atoms emerged timidly in the realm of quantum optics as 
potential photon sources, eventually alternative to their well established natural counterpart 
[T]. Nowadays, due to progress in material fabrication and characterization techniques, they 
have become not just plausible, but in cases even more suitable for both, basic research 
and technological applications. One outstanding example is the initially unambiguous and 
now impressively diaphanous observation of dressed states in self assembled quantum dots 
[2HI]. This makes nanostructured qubits (NQ) most attractive systems for doubly dressing 
endeavors [3HH]. 

On the other hand, key properties of optoelectronic systems and devices such as emission 
[absorption] frequencies and intensities, are directly determined by the local density of optical 
states (LDOS); which basically describes how efficient a system is to emit [absorb] photons 
at some particular energy. In those devices, the LDOS is set by the electronic density 
of energy states, so that reduction in the dimensionality rises sharpness in the ranges of 
optically accessible energies. These changes achieve a limit when all the spacial dimensions 
are comparable to the wavelength associated to the confined charge carriers (OD systems), 
case in which the LDOS distribution consists of Dirac deltas (actually narrow Lorentzians, 
as in high quality NQs) [9]. 

In this work, a size-independent method to control the LDOS of a OD semiconductor emit- 
ter is described; thus providing an indeed realizable scheme to generate photonic subbands 
of tunable width from fully discretized states. 

The paper is organized as follows; first a general formulation to describe slightly detuned 
double driven NQs and simulate their resonance fiuorescence spectra, is presented. Then, 
special attention is paid to the monochromatic double dressing case in which photonic sub- 
band generation is demonstrated. In the last part an analytical expression is derived for the 
relevant case under study, and conclusions are drawn. 

2. THEORETICAL DESCRIPTION 

The NQ is modeled as a semiconductor quantum dot (QD) with well defined vacuum and 
S-type exciton states, separated by energy hwo, under continuous stimulus by two lasers of 
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frequencies uja and ujb, respectively [see figure 1(a)]. 

The total Hamiltonian for the two distinguishable fields interacting with the exciton in 
the dot, in the Jaynes-Cummings framework reads 

H = -huQCTz + huA^n + 2^ %(^^+ + c-c.) 

+hwB{rh + ^) + + c.c.) + Eoff , (1) 

where (T+ is the rising part of the dot dipole transition operator, g is the dot-field coupling, 
and n and a {m and b) are the number of photons and photon annihilation operators for 
the laser A (B), respectively [IDl HI]. Eoff is an offset than allows the energy reference to 
be chosen at convenience. 

Currently, highly confined quantum dots exhibit typical neutral exciton excitation energy 
at the eV scale [12], so that detunings 6a = wq— and 6b = wq— at the order of terahertz 
or even far infrared, can be considered much smaller than the resonance frequency Uq. 

Let us assume Uq > ua > i^b and take as basis the triple direct product between the 
energy eigenstates of the dot and each of the fields, i.e. | k,n,m) where k = 1,2 and 
n,m = 0,1,2,3... It can be noted that if the coupling g is artificially turned off, for every 
fixed integer K > all states with k = and n + m = K, alongside of all states with k = 1 
and n + m = K — 1] generate a cluster with 2K + 1 eigenenergies ranging from Ei = Khcus 
(corresponding to | 0, 0, K)) to E2K+1 = {K — l)hjjA + (corresponding to | 1, — 1, 0)). 

Now, if the coupling effects are taken into account, the off diagonal non zero terms are 

(m, n, 1 I I 0, n + 1, m) = hgy/n + 1 , 

(m, n,l \ H'^ \0,n,m+ 1) = hgy/m + 1 , (2) 
and their corresponding complex conjugates. 

For a given K, we can rewrite the basis elements in the form | k, K — 171,771) , which 
emphasizes the number of photons in the laser B as compared to the total of photons. 
We introduce the label / for the elements of the basis in the corresponding cluster; i.e. 
I = 1, 2, 3, 2K, 2K + 1. Thus we can notate by | (f)f) the basis of a determined subspace 
with total number of photons K (see inset in Fig. 1). The index / is related to the number 
of photons in the laser B and to the QD state through 
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k = 0{k = l) A m = ^- ^ ; If / odd (even) . (3) 

Once the Hamiltonian in eq. (1) is written in the basis | 0f^), it can be numerically diag- 
onalized to obtain 2K + 1 eigenvalues with their corresponding eigenvectors. We represent 
with I the orthonormal basis formed by the eigenstates of the coupled system; that is 

2K+1 

I ^f) = E I <^f) ■ (4) 

The coefficients c^i = {(p^ \ ipf) , are directly obtained from the columns of the unitarian 
matrix P that diagonalizes the Hamiltonian, in such a way that P^^HP = H' (where H' is 
diagonal). They evidently represent the projection of the states of the uncoupled basis on 
the eigenvectors of the coupled system. 

When two lasers with specified powers are applied on the dot, the numbers of photons 
in lasers A and B interacting with the dot states are set; let us say and M, respectively 
(A^ + M = K). This means that the optically accessible states of the coupled system are 
those in which the specific state | 0, A^ — M,M) =\ (p^M+i)^ P^^^ '^^ superposition. 
Hence, under such conditions the LDOS of the coupled system associated to the reference 
state I is 

Pm{E) = I (<^^M+i I ^f) P S{E - En . (5) 

We are interested in the strong coupling regime, where the Rabi splitting as compared to 
the QD emission linewidth is large enough to allow steady exciton population. In this case 
the inelastic part of the light-matter scattering is the dominant one fELJ, and the emission 
intensity in a particular frequency u will be proportional to the transition rate of the system 
releasing a photon of energy hu. If the transformation P is known [and consequently the 
LDOS p^{E)], the fiuorescence spectrum can be obtained from the Fermi golden rule (FGR) 
[IMS]. Namely 

I{lo) (X ^ K^f I a_ I T/^f +1) 1^ pUEF)5{hiu + Ep - Ei) , (6) 
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where {il)p \ a- is the interaction matrix element in which the lowering operator 

turns the initial state | /) into the final state | F)\ and p^{Ef) stands for the availability 
of the final state in the I F transition. 

In order to include the finite exciton lifetime for simulating realistic spectra, the Dirac 
delta distribution in the FGR can be replaced by a Lorentzian distribution L{hu) — Ej — 
Ep;Tjp) = ~ (nuj-Ei+EFf+{hri f)^ ' where Tjp is the I F transition spectral linewidth 
(Wigner-Weisskopf approximation) [161 E] 

For an experimentally set total number of photons "i^" , the dominant emitting transitions 
are those from initial states in the cluster corresponding to the subspace + 1 to final states 
in the cluster corresponding to the subspace K. 

The relevant matrix element are then those of the form 

2K+1 2K+3 

(^f I ^- 1 = E E 4,F^(€ I ^- 1 ' (7) 

ii=i i2=i 

where the "*" in the coefficients of the final states has been omitted because the elements 
of the transformation P are all real. 

Considering the action of the operator crl, the correspondences shown in eq. ([s]), and 
the orthonormality of the basis | (f)f)', the required dipole transition matrix elements can be 
expressed in the compact form 

j=2,4,6,... 

3. MONOCHROMATIC DOUBLE DRESSING 

From this point we focus on the resonant monochromatic double dressing case, in which 

Given this condition, eq. ([T]) for each subspace becomes a tridiagonal matrix with constant 
diagonal term Khwo (a manifold with degree of degeneracy 27^^+1, if the g coupling is turned 
off). This resembles a closest-neighbor tight-binding Hamiltonian with position depending 
hopping [18], whose diagonalization straightaway renders formation of energy bands. 

Therefore, analogously to the pass from completely discretized atomic states to band 
structures in crystals; it is natural to expect a drastic change in the LDOS, evolving from 
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FIG. 1. (a) Schematics of the system and involved photon energies, (b) Energy level structure 
in different manifolds before coupling (the "hopping terms" are indicated), for the resonantly 
monochromatic double dressing case. The inset shows the top-left corner of the Hamiltonian for a 
manifold with K>3, where the corresponding basis element and their labels are indicated. 



fully discretized levels to quasi-continuous photonic subbands. 

Figure 1(b) depicts the manifolds formed due to the presence of the two lasers. The 
couplings between levels are presented in a analogous way to the hopping terms in a first- 
neighbor model for one dimensional atomic chains. 

Diagonalization of this Hamiltonian [see inset in figure 1(b)] produces 2K+1 eigenenergies 
{Ef^ = hgef) which are distributed symmetrically around 0, along the interval —hg\/2K < 
Ef^ < hgV2K. 

Figure 2(a) shows a graphical representation of the matrix transformation P which di- 
agonalizes the Hamiltonian of eq. (1) in the resonant monochromatic double dressing case, 
while figure 2(b) presents the LDOS as function of the energy and number of photons in the 
laser B (for a given K = 200). They clearly illustrate how the band structure appears as a 
direct consequence of the presence of laser B, so that the well defined side peaks (MoUow 
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FIG. 2. (a) Array-plot of the unitary transformation P than diagonahzes the Hamiltonian H. (b) 
LDOS as function of the energy and number of photons in the laser B (m). {K = 200. Energy has 
been offset, so that KHloq = 0). 

triplet pl-ffl 119]) turn into energy sidebands as the number of photons in that laser increases 
(upper and lower subbands, respect to the initially degenerated eigenenergy E = 0). 

In the purely radiative decay limit, for an exciton lifetime r, the spectral linewidths are 
respectively; Tc = 27 ^'^^ ^ transition between the upper-upper and lower-lower subbands 
(central peak), and Tg = ^ for a transition between the upper-lower and lower-upper 
subbands (side peaks) jTT| [19]. 

4. ANALYTICAL APPROXIMATION 

For the sake of physical insight and simplicity in calculations, an analytical approximation 
to obtain the transformation coefficients from well known functions can be pursued. 
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To do this, for the K-th manifold we change the basis to symmetric and antisymmetric 
combinations of the elements of the original one; i.e. | = -^{\ 0,K — m, m)± | 1,K — 

m — Because the manifold has odd number of states, the last one in the original 

basis (I 0,0, -fC)) is intentionally left unpaired for ordering the new basis in reference to it. 

In this new basis, the dot-field interaction matrix element (0^ _|_ | \ (p^, = (0^ _|_ | 
a+{a + b) + (j_(a^ + b^) \ 4^^i^±), can be further evaluated to obtain 

(0^,± I I 05,±) = ±\{26^,m'^^K^ 

+5m,m'-lVm/ + Sm,m'+lVm' + 1) 

-Sm,m'+Wm' + 1) . (9) 

If the new basis is ordered in the form 

{| 0^1,+),! 0,0,i^),| 0f _!,_), I 0f-2,-),-,| C)}- 

Therefore, the Hamiltonian takes the explicit form 
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(10) 



In the above matrix, the effects of the laser A on the QD are contained in the diagonal, 
where energy splittings of magnitude iS.E = 2hg^/n are observed between the symmetric 
and antisymmetric combinations | <j)^^) and | 

Taking aside the diagonal part , the diagonalization of the remaining matrix (H^^) 
delivers the energy modifications caused by the laser B on the dressed states of the coupled 
system QD-laser A. Those eigenvalues and eigenvectors define the energy spectrum and 
LDOS of the double driven QD, symmetrically distributed around hgy/N and —hgy/N 
(upper and lower energy subbands). 
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Matrix H can be seen as composed by four equally sized blocks, plus the row and 
column corresponding to the state | 0,0,7^). The off-diagonal terms in the top-right and 
bottom- left blocks (off-diagonal blocks, responsible of the mixing between the diagonal ones), 
increase in absolute value as their position are closer to the matrix center, i.e. m ^ K. 

Hence, for the system in the limit K ^ oo and M / — )■ (in which laser A is much more 
intense than laser B and then mixing effects between blocks are negligible); the eigenvectors 
from diagonalization of H^^ are a superposition of the new basis according to 

\^^.±) = fl~4i,±\~^l±) ' (11) 

i=i 

where the coefficients cfi ± = {(pf^± \ 'ipi^±) are to be associated to a set of orthonormal 
functions. 

Since the off-diagonal elements in the top-left [bottom-right] block, located exclusively 
right above and below the matrix diagonal, have the same structure as those of the operator 
(b + 6^) = y/2x [—(b + b'^) = —\^x] written in the eigenbasis of the number of photons 
operator (m = b'^b) [2D]; in this limit the matrix elements of the transformation that diago- 
nalizes H^^ can be obtained in good approximation from the harmonic quantum oscillator 
eigenfunctions $j (z) [ZllIIlEO]. This is 

cf,i,±^cj,e,± = {j\±e)^^,i^) , (12) 

with e = \/2x (— oo < e < oo), a dimensionless continuum parameter which multiplied 
by hg becomes the eigenenergy corresponding to the eigenvector | \/2x). The functions are 
explicitly 



where Hj{x) is the j-th order Hermite polynomial [121 1^ . 

For the relevant transitions between each of the two subbands in the K + 1 and K 
manifolds, there are four possibilities: transitions between the upper (lower) subbands, 
alongside with transitions from the upper K + 1 (lower -|- 1) subband to the lower K 
(upper K) one. The dipole transition matrix elements become | {ipp± \ \ P= 
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±|5(£_F — ei) and | {'4'f,t I ^- I 1^^ ^^H^i + ± 2-\/iV), respectively; where the 

transition energy has been also shifted to make hwQ = 0. 

On the other side, the normalized LDOS for the reference number of photons M as 
function of the energy of the transition final state {Ep = hgep), is now given by 



Pm(E,) = (^) I 4«(^, 



(14) 



Figure 3 shows the LDOS as function of the emission frequency normalized to the coupling 
g, for different numbers of photons in the laser B. It is worth noting how the LDOS has 
absolute maxima at values close to a; = ±v^M(yf, respectively [21]. For frequencies more 
separated from the reference {<^/g = 0), it decays rapidly to zero; setting an approximate 
width for the subbands of two-times the Rabi splitting associated to the laser B. 




FIG. 3. LDOS for different numbers of photons in the laser B. The thin vertical lines indicate the 
corresponding values oj = ■^\fMg. The dashed black line shows the case in which only laser A is 
applied (M = 0). 

Under these considerations, the FGR yields for the emission spectrum 



OO /"OO 



/(w) oc / \ dEi dEp \b{Ei + Ep + 2hgVN) + 26{Ei - Ep) 

J —OO J —OO 

+S{Ei + Ep- 2hgy/N)] puiEp) L{huj + Ep - Ej, Tj^p) 



(15) 



which after insertion of eq. ( 14 ) and integration over Ej turns in 
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FIG. 4. Calculated resonance fluorescence spectra for different laser intensity ratios; given a fixed 
K = 100. The solid lines are obtained from eq. ^ while the dashed lines are obtained from eq. 



(16). The dotted black line shows the well known Mollow triplet case, when just laser A is applied. 



(a) Tc = 0.5(7 s-iicl = 0.75^, (b) Fc = 5g and Fs = 7.5g. 

In figure 4 the resonance fluorescence spectra calculated for different ratios between the 
numbers of photons in lasers A and B (a^ = are shown for a fixed total number of 
photons K = 100. They are obtained by using both, the numerical quasi-exact approach 
of eq. ^ and the analytical approximation of eq. (16). In the upper frame, for a close 



to the value of the coupling constant g (i.e. the exciton lifetime is large enough respect to 
the Rabi oscillation periods defined by the laser intensities); it is clearly appreciated how 
the Mollow side peaks of the single excitation laser case evolve into sidebands, reflecting the 
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energy band formation under the double dot-field coupling. Markedly the subband width 
increases with the number of photons in the laser B. 

In other words, due to the second driving field, well defined peaks spread into optically ac- 
tive regions; then exhibiting similarities with LDOS proper of higher dimensionality systems 



When Tc is substantially larger than g, the effects of the second laser become less no- 
ticeable and useful. This because in the time dominion, the short exciton lifetime inhibits 
coherent Rabi oscillations making almost irrelevant the presence of laser B. This sensitiv- 
ity of the system to the ratio between g and Fc provides a way to estimate the order of 
magnitude of the dot-field coupling. 



Figures 4(a) and 4(b), evidence that the analytical approximation behind eq. (16), as 
long as the value of a is not close to one, works ostensibly well even for a moderate number 
of photons. The main discrepancy is found around the subband edges, where the exact 
calculation predicts slight asymmetry due to the finite nature of the Gilbert subspaces. 
Whereas such an asymmetry is tenuous, actually it has been experimentally observed [8]. 



5. CONCLUSIONS 

In conclusion, a theoretical approach to model and simulate doubly driven artificial atoms 
was implemented. 

As a main result, tailored manipulation of the optical density of states in semiconductor 
quantum dots, has been shown for the case of monochromatic double dressing. It was 
described how by coupling a nanostructured qubit simultaneously to two distinguishable 
lasers whose frequencies match the exciton transition, a discrete eigenstate turns into an 
energy subband in a process closely analogous to band formation in solid state physics. 

Such strong changes in the local density of optical states, controllable through the ratio 
between the laser intensities; open new possibilities for on-demand photon emission from 
artificial atoms. 

The presented results are in remarkable qualitative and quantitative agreement with 
experimental measurements, as presented in Ref. [8]. 
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